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ABSTRACT 


This  report  presents  a  method  to  design  the  inviscid  wall  contour 
of  a  supersonic,  axisymmetric  nozzle  producing  uniform  parallel  flow 
at  the  exit.  The  calculation  methods  used  are  restricted  to  a  perfect 
gas.  The  axis  velocity  distribution  must  be  specified.  A  complete 
flow  pattern  within  the  nozzle  is  computed  by  the  method  of  character¬ 
istics  beginning  at  Mach  1. 0  on  the  axis  and  extending  to  the  region  of 
uniform  parallel  flow  at  the  exit.  The  wall  points  are  calculated  by 
means  of  a  mass  integration  technique.  The  methods  of  calculation 
included  have  been  programmed  in  Fortran  II  language,  and  a  complete 
listing  of  the  program  is  given  in  Appendix  I. 
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NOMENCLATURE 


Distance  along  a  characteristic 
Velocity 

Velocity  /limiting  velocity  (V/Vmax) 
Physical  abscissa 
Physical  ordinate 

Mach  angle  (sin  a  -  1. 0/Mach  number) 
Ratio  of  specific  heats 
Inflection  angle 

Flow  anple  (with  respect  to  axis) 
Density 

Prandtl-Meyer  angle 


SUBSCRIPTS 


Refers  to  right-running  characteristic 
Refers  to  left-running  characteristic 
Origin  of  right-running  characteristic 
Origin  of  left -running  characteristic 
Computed  point 

End  points  on  line  BC 

Normal  Component 
Point  on  line  BC 
Point  on  line  CD 
Total  (stagnation  condition) 
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1.0  INTRODUCTION 


This  report  presents  a  method  to  design  the  inviscid  wall  contour 
of  a  supersonic,  axisymmetric  nozzle  producing  uniform  parallel  flow 
at  the  exit.  The  calculation  methods  used  are  restricted  to  a  perfect 
gas.  The  axis  velocity  distribution  must  be  specified. 

A  complete  flow  pattern  within  the  nozzle  is  computed  by  the  method 
of  characteristics  beginning  at  Mach  1. 0  on  the  axis  and  extending  to 
the  region  of  uniform  parallel  flow  at  the  exit.  The  wall  points  are 
calculated  by  means  of  a  mass  integration  technique. 

The  methods  of  calculation  included  in  this  report  have  been 
programmed  in  Fortran  II  language,  and  a  complete  listing  of  the  pro¬ 
gram  is  given  in  Appendix  I.  This  program  has  been  run  on  the 
IBM  7074  computer  at  the  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Systems  Command  (AFSC),  to  obtain  the  results 
to  be  presented. 

2.0  SCOPE  OF  CALCULATIONS 


For  the  purpose  of  calculations  the  flow  in  the  nozzle  is  arbitrarily 
divided  into  three  distinct  regions  as  shown  in  Fig.  1 . 


Origin 

of 

Source 


Fig.  1  Major  Nosxla  Divisions 
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2.1  THROAT  REGION 


In  calculating  flow  within  the  throat  region,  points  A  and  B  of 
Fig.  2  are  given  points. 

The  velocity  distribution  on  AB  is  given  either  by  an  equation  or  as 
a  series  of  discrete  points.  The  flow  angle  at  point  C,  the  inflection 
point,  is  also  an  input  value.  Knowing  this  angle  and  the  velocity  at 
point  B,  a  left -running  characteristic  BC  (the  upstream  boundary  of 
the  source  flow  region)  can  be  established  using  the  source  flow 
relation 


0p  =  +  20P 

where  the  subscript  p  refers  to  a  point  on  the  characteristic  line  BC. 
Line  BC  is  then  divided  into  an  arbitrary  number  of  segments  having 
an  equal  A0  between  successive  points.  These  are  initial  points  for 
right-running  characteristics  from  line  BC.  The  mass  flow  in  the 
nozzle  is  found  by  integrating  along  BC  using  the  technique  outlined  in 
Section  2.4. 

Ar  arbitrary  number  of  points  are  specified  as  the  origin  points 
for  the  left -running  characteristics  beginning  on  the  axi  a.  Character¬ 
istic  lines  are  calculated  until  a  field  point  is  computed  from  the  last 
known  wall  point.  As  the  last  characteristic  point  lies  beyond  the  wall, 
a  mass  integration  technique  is  applied  to  find  the  wall  point  location 
on  the  characteristic  line. 

22  SOURCE  FLOW  REGION 

The  area  enclosed  by  BCD  in  Fig.  1  consists  entirely  of  source 
flow,  and  no  characteristic  solutions  are  computed. 
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22  DOftXSTRSAM  REGION 


Line  CD  is  a  right  -  running  characteristic  line  on  the  downstream 
boundary  of  the  source  flow  region.  This  characteristic  is  calculated 
in  a  manner  similar  to  line  BC  in  the  throat  region.  Here  the  relation 
used  to  establish  this  line  is 


-  lAc  +  2(j?  ~  V 


In  a  manner  analogous  to  that  used  in  defining  the  axis  line  AB,  it 
is  necessary  to  designate  points  D  and  E  and  to  specif'  an  axial  velocity 
distribution  to  establish  line  DE.  Also  similar  to  the  pi  ocedure  used 
for  Region  I,  line  DE  is  divided  into  an  arbitrary  number  of  segments. 
The  axis  points  are  the  origins  of  right- running  characteristics. 

Successive  right-running  characteristics  originating  from  line  DE 
are  computed,  and  wall  points  are  found  by  mass  integration  as  in  the 
throat  region. 

Line  EF  in  Fig.  3  is  a  left- running  characteristic  consisting  of 
uniform  parallel  flow  at  the  desired  design  exit  Mach  number.  Line  EF 
is  a  straight  line  divided  into  an  arbitrary  number  of  equal  intervals  for 
origin  points  to  right -running  characteristics.  The  mass  integration 
along  characteristics  from  line  EF  is  handled  in  a  different  manner 
than  integrations  beginning  on  the  axis.  Here  the  calculations  are 
performed  with  an  initial  mass  integration  value  of 

mo  15  (7—-)  X  Total  Maas 
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2A  MASS  INTEGRATION  TECHNIQUE 

The  differential  equation  governing  mass  flow  at  a  point  on  a 
characteristic  line  is 


dM  -  2  it  •  r  •  Vfl  •  p  •  ds 


(1) 


The  normal  component  of  velocity  at  a  point  on  either  a  right -running 
or  a  left-running  characteristic  is 

Vjj  ■  V  •  sin  a 


a.  Left -Running 


Fig.  4  Characteristics 


y-i 


and  V 


Let  p  -  {^)  ’  p'  where  =  (1  -  W*) 

E  •  VmM.  Now  Eq.  (1)  may  be  expressed  as 


v... 


r  BIX 


y  _  l 

dM  *  2w  •  r  •  W  •  Vmmx  •  sin  a  •  (1  -  W*)  •  px  •  ds 


(2) 


For  a  given  nozzle  condition,  the  factor  2^  •  VBIX  •  pt  is  constant.  By 


letting 


•  V..t  •  pt 

nondime  ns  ionalized  mass  differential  equation 


and  substituting  in  Eq.  (2),  one  obtains  the 


dm  -  r  •  f  •  (1  -  W*)y“1  . 


sin  a  • 


ds 


(3) 


By  defining 


r  •  W  •  (1  -  W1) 


y-i 


sin  a  *  f  (s) 
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Eq.  (3)  can  be  simplified  to 


dm  •  f(a)  ds  (4) 

To  calculate  the  nondimensional  flow  in  the  nozzle,  the  following 
technique  is  used. 

In  the  throat  region,  the  left -running  characteristic  defining  the 
origin  of  the  source  flow  region.  Region  II  in  Fig.  1,  is  known  from  the 
axis  to  the  wall.  Equation  (4)  is  integrated  as  follows: 

For  each  point  on  the  line  BC,  the  value  f(s)  is  calculated.  The 
distance  (s)  from  the  axis  to  any  point  is  the  sum  of  the  straight  line 
segments  from  point  to  point. 

The  mass  integral  is  computed  by  analytically  evaluating  a  para¬ 
bolic  fit  through  each  three  consecutive  points  and  summing  the  integral 
values.  The  limit  of  integration  of  each  integral  is  from  the  first  to 
the  second  point. 

Computations  for  mass  in  the  downstream  region.  Region  III  in 
Fig.  1,  is  carried  out  in  a  similar  form  except  that  integration  is  along 
the  right-running  characteristic  line  CD. 

2.5  LOCATING  WALL  POINTS  ON  A  CHARACTERISTIC 

The  final  point  calculated  on  each  characteristic  line  will  be  beyond 
the  corresponding  wall  point  for  that  line.  Hence,  if  the  nondi  men  - 
sionalized  mass  flow  is  integrated  to  the  final  point  it  will  exceed  the 
nozzle  mass  flow.  To  find  the  unknown  wall  point,  the  sum  of  the  mass 
flow  from  point  to  point  is  calculated  until  a  point  on  the  characteristic 
line  is  reached  where  the  line  mass  integral  exceeds  the  known  nozzle 
mass  flow.  An  iteration  technique  is  applied  to  find  the  unknown  wall 
point.  The  values  of  the  wall  point  are  calculated  by  linear  interpolation 
of  the  known  end  point  values. 


3.0  RESULTS 


A  comparison  of  theoretical  results  and  the  characteristic  solution 
in  the  transonic  region  is  shown  in  Fig.  5.  The  points  shown  were 
calculated  by  the  program  included  in  this  report. 
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One  important  factor  in  the  proper  use  of  this  program  is  the 
selection  of  the  initial  points  of  the  characteristic  line  beginning  on  the 
axis.  The  spacing  will  depend  upon  the  axis  velocity  distribution.  It 
has  been  found  that  a  value  of  1. 5  <  p  <  2.0  will  result  in  a  good  spacing 
in  the  throat  where 

Xt  =  Abscissa  at  Mach  number  1.0 

X,  =  Abscissa  at  the  beginning  of  source  flow  region 

N  =  Nth  point  1  £  N  £  M 

M  =  Total  number  of  points  on  the  throat  axis 
X  =  Abscissa  corresponding  to  the  Nth  point 

X  -  Xt  =  /N  -  1  \p 
X,  -  X,  \M  -  1/ 

A  prior  computer  program  was  written  to  obtain  the  inviscid  wall 
contours  by  means  of  an  extrapolation  technique  based  on  the  starting 
known  wall  point  and  extending  this  along  the  wall  streamline  by 
extrapolation  of  the  flow  direction.  The  technique  was  found  to  be 
unsatisfactory  since  errors  tended  to  be  propagated  from  point  to  point. 
It  was  a'.so  found  that  the  network  size  greatly  affected  the  accuracy  in 
locating  the  wall.  With  the  mass  flow  integration  technique  used,  no 
propagated  error  occurs  since  the  calculation  of  each  wall  point  is 
based  on  a  single  characteristic  line.  The  total  number  of  points  needed 
in  the  characteristic  field  to  obtain  satisfactory  inviscid  wall  contours 
was  considerably  less  than  required  with  the  streamline  extrapolation 
technique. 


4.0  CONCLUSION 


A  comparison  was  made  of  the  flow  field  in  the  throat  region 
calculated  by  the  techniques  in  this  report  to  analytical  results  of  other 
sources.  Excellent  agreement  of  results  was  found  between  the  methods 
(see  Fig.  5). 

With  this  program  an  inviscid  wall  contour  can  be  designed  with  an 
initial  wall  point  on  the  left-running  characteristic  corresponding  to 
Mach  1.0  on  the  axis  to  an  exit  having  uniform  parallel  flow.  It  must 
be  kept  in  mind  that  the  methods  applied  hold  for  a  perfect  gas  and  that 
the  axis  velocity  distribution  must  be  specified. 
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O  Characteristic  Points  from  Program 

A  V/V  from  Program 
max 

—  Ref.  2 

—  Ref.  3 


0  0.05  0.10  0.15 
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Fip.  5  Comparison  of  Thooroticol  o nrf  Computed  Results 
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APPENDIX  I 

COMPUTER  PROGRAM  FOR  NOZZLE  CONTOUR 


The  program  presented  here  is  based  on  the  methods  outlined  in 
the  preceding  sections.  As  previously  pointed  out,  there  is  a  difference 
in  the  manner  of  computing  the  throat  and  downstream  regions,  and 
they  may  be  computed  cither  separately  or  together. 

In  either  case,  a  complete  set  of  input  data  must  be  prepared  for 
the  region  under  consideration.  However,  if  the  downstream  region 
is  computed  immediately  after  the  throat  region,  the  X-location  of  the 
inflection  point  need  not  be  specified  as  it  is  calculated  as  the  wall 
point  on  line  BC  in  Fig.  1. 

Care  must  be  exercised  in  specifying  the  distance  between  points 
on  the  axis  to  ensure  that  the  characteristic  points  to  be  computed 
will  be  properly  spaced.  If  axis  points  A1  and  B1  in  Fig.  6  are  chosen 
judiciously,  A2  will  be  computed  in  a  manner  to  make  it  almost 
equidistant  from  its  parent  points.  This  effect  will  be  propagated 
throughout  the  field. 
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Fig.  6  Axis  Distribution 

The  number  of  points  on  an  initial  characteristic  line  in  either 
region  should  be  chosen  such  that  the  dimension  size  in  the  program 
is  not  exceeded  as  points  are  added  to  successive  characteristic  lines. 
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NOMENCLATURE  OF  PROGRAM 


A 

B 

C 

S 

FS 

AXIS 

BEGIN 

FINAL 

WALL 

IEND 

M 

N 

IX 


IWOT 


IP 


NP 


Parameters  associated  with  right-running  characteristics 
Parameters  associated  with  left- running  characteristics 
Coefficients  of  polynomial  describing  the  axis 
Characteristic  arc  length 

Integrand  value  used  in  mass  integration  technique 

Refers  to  parameters  located  on  the  centerline 

Refers  to  parameters  on  the  initially  calculated  character¬ 
istic,  i.  e. ,  line  BC  in  the  throat  region  and  CD  in  the 
downstream  region. 

Refers  to  the  parameters  calculated  at  points  along  line  EF. 

Refers  to  the  values  calculated  at  the  wall  points  found  by 
the  method  described  in  Section  2.5. 

Control  to  end  program  or  to  calculate 

Number  of  points  taken  on  BEGIN  line,  BC  or  CD 

Number  of  axis  points 

Control  dependent  on  method  used  to  set  up  axis 

IX  =  -1  if  discrete  X-coordinates  are  to  be  read  in  preceded 
by  the  corresponding  values  of  V/Vmax;  coefficients, 
C,  not  used. 

IX  =  0  if  discrete  X-coordinates  are  to  be  read  in  and 

V/Vmax  is  to  be  calculated  from  the  coefficients,  C. 

IX  =  +1  if  axis  points  are  to  be  computed.  Coefficients  are 
to  be  read  in. 

Control  on  amount  of  calculated  points  to  be  written  out. 

IWOT  =  0  writes  out  wall,  axis,  and  first  and  last 
characteristics 

IWOT  i  0  writes  out  all  calculated  points 

Indicates  which  region  is  to  be  computed 

IP  =  0  for  Region  I 

IP  i  0  for  Region  III 

Number  of  points  to  be  taken  on  line  EF 
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GAM  Ratio  of  specific  heats 

ETA  Angle  at  inflection  point,  point  C,  in  degrees 
XC  X-coordinate  at  point  C 

XI  &  X2  V/Vmax  is  a  polynomial  dependent  on  an  Xr  relative  to 
source  How  such  that 


v/vmax 


w  =  «5o 


where 


Xr  -  XI 
X2  -  XI 


P  Governs  distribution  on  axis.  P  =  1  for  linear  distribution. 

P  >  1  packs  axis  points  toward  low  end  of  axis. 

C  Constants  used  in  defining  the  velocity  ratio  along  the  axis 

in  terms  of  V/Vmax. 

XHI  Maximum  X-coordinate  on  the  axis,  point  B  or  E 

XLO  Minimum  X-coordinate  on  the  axis,  point  A  or  D 

AXIS(1,  N)  Axis  coordinates  to  be  read  in  when  IX  ^  0 

AXIS(3,  N)  V/Vmax  corresponding  to  axis  points  to  be  read  in  when 
IX  =  -1 


Subscripts  (on  A,  B,  AXIS,  BEGIN,  FINAL,  and  WALL) 

First  subscript  (1-6)  1  =  X-coordinate 

2  =  Y- coordinate 

3  *  V/Vmax 

4  =  Mach  angle 

5  =  Flow  angle 

6  =  Mach  number 

Second  subscript  (1-100)  designates  point  number. 
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Initial  Block 
Dimension  statement 
Formats 

Alpha  characters 
Test  to  compute 


or  stop 


400  Block 


Calculate  points  on 
first  characteristic 
taken,  both  regions 
Compute  mass 


500  Block 

Calculates  charac¬ 
teristics  origi¬ 
nating  on  the 
axis 


Writes  out  vail 

points  calculated 
in  700  Block 
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INPUT  DATA  AND  SAMPLE  VALUES 

(Input  values  of  IEND,  M,  N,  IX,  IWOT,  IP,  NP.  GAM,  ETA,  XA, 
XI,  X?,  P,  C,  XHI,  XLD,  and  AXIS  discussed  in  Nomenclature  of 
Program. ) 

Sample  input  used  for  design  criteria  of  Mach  8,  nozzle  half  angle, 
ETA  =12  deg. 

Throat  region  (Region  1): 

IEND  =  1  (to  proceed  with  calculations) 

M  =  15  (15  right-running  characteristics  from  line  BC) 

N  =  20  (20  axis  points  right-running  characteristics) 

IX  =  1  (to  calculate  distribution  along  axis) 

IWOT  =  1  (to  write  out  all  characteristics) 

IP  =  0  (to  calculate  in  the  throat  region) 

NP  =  0  (to  calculate  in  the  throat  region) 

GAM  =1.4 
ETA  =  12  deg 

XC  =  0.0  (not  used  for  throat  calculations) 

XI  =  0.4727  (discussed  in  Nomenclature  of  Program) 

X2  =  1.91607  (discussed  in  Nomenclature  of  Program) 

P  =  2.0  (square  distribution  on  axis) 

C(l)  =  0.408248247 
C(2)  =  0.551273121 
C(3)  =  0.010438091 
C(4)  =  -0.369697773 
C(5)  =  0.229213428 
C(6)  =  -0.042703583 
XHI  =  1.91607 
XLO  =  0.4727 

Downstream  Region  (Region  HI): 

IEND  =  1 
M  =  40 
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N  =  6 
IX  =  1 
IWOT  =  1 
IP  =  1 
NP  =  35 
GAM  =  1.4 
ETA  =  12  deg 

XC  =0.0  (calculated  in  throat  region  and  not  necessary  to  read  in) 

XI  =  13.015491 

X2  =  15.211972 

P  =  1.0 

C(l)  =  0.961279408 

C(2)  =  0.005421945 

C(3)  =  0.005421945 

C(4)  =  0.001807315 

C(5)  =0.0 

C(6)  =0.0 

XHI  =  15.211972 

XLO  =  13.015491 

IEND  =  0  (to  end  program). 
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05314  MAIN  PROGRAM  AXI-SYMMETR IC  SUPERSONIC  N022LE 
PROGRAMMED  BY  WC  MOGER  AND  DB  RAMSAY 

COMPUTATION  BASED  ON  MASS  FLOW  INTEGRATION  BY  PARABOLIC  FIT 


DIMENSION  A<6»100) *B(6»l00>  .  C<6)  *D< 10) *FS( 100) *S( 100) 
DIMENSION  AX  I S<6*l00) .BEGIN (6*100 ).FINAL< 6*100 ) *WALL( 6*100) 
EQUIVALENCE  ( BEGIN! 1 ) » FINAL ( 1 ) ) 

1  FORMAT  (9H  FN  05314) 

2  0  FORMAT  (I1.78H 

2  1  ) 

3  FORMAT  (1615) 

4  FORMAT  (6E12.0) 

5  FORMAT  (/////10X.5HINPUT///5X* 

1  45HNUMBER  OF  POINTS  ON  FIRST  CHARACTERISTIC  (M)«I3// 

2  5X.29HNUMBER  OF  POINTS  ON  AXIS  (N)*I3//5X.6HGAMMA*F7.4// 

3  5X • 2 3H INFLECT  ION  ANGLE  ( ETA ) -F7.3 *2X .7HDEGREES// 

4  5X*36HCOORDINATE  OF  INFLECTION  POINT  (XC)«F9.5// 

5  5X.38HFACTORS  IN  V/VMAX  TO  DETERMINE  X  ( XI ) -F9.5//38X 

6  5H(X2)*F9.5//5X.38HPOWER  GOVERNING  AXIS  DISTRIBUTION  <P)»F8.5 

7  //5X» 31HCOEFF IC IENTS  IN  TERMS  OF  V/VMAX// 10X .5HC ( 1 ) «1PE13. 7 

8  //10X#5HC(2)*lPE13*7//10X.5HC(3)=lPE13.7//l0Xt5HC(4)«lPE13.7 

9  //10X.5HC(5)«1PE13.7//10X#5HC(6)=1PE13.7) 

7  0  FORMAT  ( 1H1 *10X *A5// 10X .5HPOINT f8X . 1HX. 14X. 1HY . 12X *6HV/V«AX 
7  1  6X  1 1 2HMACH  ANG. ( D ) *3X • 12HFLOW  ANG. ( D ) *4X8HMACH  NO.//lO<lOX 

7  2  .I3.2X.1P6E15.7/)/) 

8  0  FORMAT(lHllOXA5*I3//10X.5HPOINT.8X.lHX*l4X.lHY.12Xt6HV/VMAX 
1  6X.12HMACH  ANG. ( D ) *3X . 12HFLOW  ANG. (D) .4X8HMACH  NO.//) 

9  FORMAT  ( 10X* 1 3  »2X 1 1P6E1 5 .7  ) 

10  FORMAT  (1CHOFELD  LINE. 13. 6H  POINT*I3) 

11  FORMAT  ( 10X.I3*2X.1P6E15.7*2X.6H(MASS) ) 


15  TYPE  1 

WRITE  OUTPUT  TAPE  24*1 
XXC-1.0 

IF  I RS  «  6669798283 
IAXIS  «  61876982 
ICHAR  «  63686179 
IWALL  *  86617373 
CONV  «  57.29578 
17  READ  2 » IEND 

IF( IEND)100.804.100 

100  DO  101  J*1 *600 
AXIS(J)*0.0 

101  A(J)*0.0 
SENSE  LIGHT  0 

WRITE  OUTPUT  TAPE  24*2*IEND 
READ  3*M»N»IX  » I  WOT  * IP.NP 
READ  4*GAM*ETA»XC*X1»X2.P 
REAO  4* C 

WRITE  OUTPUT  TAPE  24*5 *MtN*GAM»ETA*XC*Xl » X2.P* I  Cl J ) . J«l*6) 

eta«eta/conv 

IF ( XC ) 200  *103  *200 
103  XC-XXC 
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200  SERIES  SETS  UP  FUNCTIONS  OF  GAMMA  USED  THROUGHOUT 
THE  PROGRAM 

200  G9  «  2*0/ ( GAM-1 *0) 

G8  •  1.0/G9 
G7  «  GAM+1.0 
G6  »  GAM-1.0 
G5  ■  G6/G7 
G4  »  SORTF ( G5 ) 

G2  ■  1.0/G4 

G1  •  1.0/G6 

EXP  •  0.5/G5 

CONST  •  ( 0.5*67 )**EXP 

EXP  «  -EXP 

NN*1 

300  SERIES  SETS  UP  DISTRIBUTION  ON  THE  AXIS 

300  IF( IX) 310.305 *301 

301  READ  4.XHI.XL0 
DX-XHI-XLO 

FN  »  FLOATF(N-l) 

DO  304  J*1  *N 
IF( IP >303*302. 303 

302  K*J 

GO  TO  304 

303  K-N-J+l 

304  AXIS ( 1 *K ) «DX* ( FLOATF (N-J )  /FN ) **P*XLO 
GO  TO  306 

310  READ  4*<AXIS<3*J> »J*1.N) 

305  READ  4*  (  AX IS( 1  * J) * J«1 *N ) 

306  DO  307  J-l.N 
X-AXIS(l.J) 

X»(X-X1 ) / ( X2-X1 ) 

IF ( IX ) 31 1 f 312  *312 

311  W»AXIS( 3  ) 

GO  TO  313 

312  W«C<1  M-X*<C<2M’X*<C<3)+X*<C<4H-X*<C<5H-X*C<6))n ) 

AX  I S ( 3.J)«W 

313  WW  ■  W*W 

XM  »  SQRTF(G9*WW/( 1.0-WWj) 

AX  IS ( 6* J )  •  XM 

S< J)»1.0/XM 

AXIS( 4* J ) »ASINF ( S( J ) ) *CONV 

307  CONTINUE 

WRITE  OUTPUT  TAPE  24* 7* IAX IS* <K *  I  AXIS! J *K ) » J»1 *6 ) *K»1 *N > 
CM*AXIS(6*N) 

DO  308  J*1 *N 

308  AXIS(4»J)»S(J) 

IF( IP  1309*400* 309 

309  XM*COSRF<ETA)/XC 
XM*XM*XM 
XM«FMA(XM*GAM) 

AXI5(6)»XM 
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400  SERIES  SETS  UP  FIRST  CHARACTERISTIC  AND 
COMPUTES  MASS 

400  PO«G2*ATANF<  S<5RTF<G5<M  AX  IS  <  6  >*AXIS<6)-1 .0))  ) 

1  *1 •5707963+AS INF  <1.0/AXIS<6)  ) 

EM  -  ETA/FLOATF  <  H— 1 ) 

DO  401  J« 1 *M 
T  «  FLOAT F ( J-l )*EM 
I F ( IP)  403*  402  *  403 

402  XM«FMV  (  PO+T-f  T  tGAM  ) 

GO  TO  404 

403  T-ETA-T 

409  XM«FMV<PO+2.0*FL0ATF< J-l )*EM.GAM) 

404  W«XM*XM*G8 

R  «  SORTF< 1.0/<CONST*XM*< < 1,0+G8*XM*XM)**EXP)) ) 

I F ( IP >406*405 *406 

405  K«J 

GO  TO  407 

406  K-M-J+l 

407  A< 1 ) «R*COSRF ( T ) 

AI2**0»R*SINRF(T) 

A(3fK)«S0RTF(W/( l.U+W) ) 

A(4»K )«1*0/XM 
AI5tM-T 

401  A (6.K ) «XM 

410  I END  «  6*M 

DO  411  J*1  * IEND 

411  BEGIN! J)  «  A <  J ) 

BEGIN ( 2*1) *0.0 
BEGI N<  5  # 1 ) »0.0 

420  CO  421  J»1 #6 

421  WALL ( J )  -  A  <  J  »M ) 

S<1)«0.0 

430  DO  431  J«2tM 
W  •  A  <  3  « J  ) 

DX«A(1*J)-A(1*J-1) 

DY«A<2tJ)-A(2»J-l ) 

S<  JI-S<  J-l  J-fSQRTF <DX*DX+DY*DY) 

431  FS< J)»A<2#J)*W*A<4# J)*<1.0-W*W)**G1 

440  LAST  -  M  -  2 
SUM  ■  0.0 

DO  441  J*1 .LAST 
K»J 

CALL  PARAB(S(K) »FS(K) tD(l ) ) 

441  SUM«P0LNT<D<1 ) *S ( K ) .S<IU1 ) )+SUM 
SUM»P0LNT(D(1  )  *5<LAST  +  1J  »S( LAST 4-2  )  )  +  SUM 
XMASS  ■  SUM 

460  DO  4S1  J*1 *M 
BEGIH<4.J)»ASINF<BEGIN<4*J) )*CONV 

461  BEGIN ( 5»J )  ■  BEGIN ( 5 * J >*CONV 

WRITE  OUTPUT  TAPE  24»7»IFIRS  *<  K  ♦<  BEGIN  <  J  )  »J«  1 .6  )  »IC«1  »M  J 
470  LAST  «  M 
LINE  ■  2 
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500  SERIES  COMPUTES  CHARACTERISTICS  ORIGINATING 
ON  THE  AXIS 

500  DO  501  J«1.6 

501  B(J)  *  AXIS(J*LINF) 

502  DO  5C6  J*1 *LAST 
X  «  J 

IF(IP)504.503.504 

503  CALL  OFELDIAI l.X) »BI l.X) *B< l.X+1) .G9) 

GO  TO  505 

504  CALL  OFELD(B( 1*X)*A(1»X)*B(1*X+1) *G9) 

505  B(6.X+1)*1.0/B(4.X+1 ) 

IF (SENSE  LIGHT  1)507.506 

506  CONTINUE 
GO  TO  508 

507  TYPE  10. LINE. J 
LAST  *  J 
SENSE  LIGHT  1 

508  LASTP-LAST+1 
IGO*l 
ME-LINE 

I F I  IP) 509. 700. 509 

509  FINAL (6. NN ) *0.0 
GO  TO  700 

538  IF (N-LINE 1539.550.539 

539  LINE«LINE+1 

DO  540  1*1 .600 

540  A( I )*B< I ) 

GO  TO  500 

541  DO  542  J*1 .LINE 

WALL ( 4. J )  *  AS  I NF ( WALL ( 4 . J ) )*CONV 

542  WALL ( 5. J ) *WALL ( 5 . J )*CONV 
I*6*LINE 

DO  543  JJ*1.I 

543  B(JJ)*WALL(JJ) 

JJ-LINE 

DO  545  L-l.LINE 
DO  544  I*!. 6 

544  WALL (I.L)*B(I»JJ) 

545  JJ*JJ-1 

WRITE  OUTPUT  TAPE  24.7. I WALL *  ( X  »  ( WALLi J.K).J*1.6).K*1*LINE) 
XXC*WALL( l.X) 

GO  TO  800 

550  IF( IP) 60 0.54 1*600 

600  SERIES  COMPUTES  CHARACTERISTICS  ORIGINATING  ON 
THE  FINAL  LINE  OF  THE  DOWNSTREAM  SEGMENT 

600  IF(NN~l)603»60l»603 

601  NNP*N+1 
W«AXIS(3*N) 

FS(NP)«SQRTF(2.0*XMASS/(W*( 1.0-W*W)**G1 ) ) 

DELX.FSI NP ) *SORTF ( CM*CM-1 .0  ) 

FN*1 • 0 /FLOAT F ( NP~1 ) 

DO  602  JJ-l.NP 
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F=FL0ATF( JJ-1 ) *FN 

FINAL  (1»JJ)*AXIS(1.N) +F*DELX 

FINAL<2.JJ)«F«FS(NP) 

FINAL! 3»JJ) «W 
FINAL (4»JJ)*AXIS(4»N) 

FINAL ( 5 . JJ ) *0.0 

602  FINAL(6»JJ)*XMASS*F*F 

603  NN=NN+1 

DO  604  1=1 #600 

604  Am=B(I) 

DO  606  J*1 *600 

606  B(  J ) “FINAL ( J  ) 

DO  607  J*NN.LAST 
K«J 

CALL  OFELD(B(l»K)»A(l»K) »Bt  1 .K  +  l )  .69) 

B(6»K+1)*1«0/B(4»K+1 ) 

IF(SENSE  LIGHT  1)608.607 

607  CONTINUE 
GO  TO  609 

608  TYPE  10.NN.J 
LAST  * J 

609  LASTP«LAST+1 
I GO=NN 
ME=NNP 

GO  TO  700 

640  NNP-NNP+1 
IF(NN-NP)600. 641*600 

641  NL=NNP~1 

DO  642  J=1.NL 

WALL(4.J)=ASlNF( WALL ( 4 . J ) ) *CONV 

642  WALL ( 5 • J ) “WALL ( 5  » J ) *CONV 

WRITE  OUTPUT  TAPE  24 . 7. 1  WALL .  MC  .  <  WALL  <  J.IC  )  »  J«1 .6)  *K«1  .NL ) 
GO  TO  800 

700  SERIES  WRITES  OUT  CHARACTERISTICS  AND  INTEGRATES 
ALONG  CHARACTERISTIC  FOR  WALL  POINT 

700  IFUWOTJ701, 703*701 

701  WRITE  OUTPUT  TAPE  24 . 8 . ICHAR .ME 
DO  702  J* I GO.LASTP 

Z  =  ASINF(BU.J)  )*CONV 

T  «  B<5»J)*CONV 

IF  I IGO-J )723»721»723 

721  IF( IP)722.723t722 

722  IF(N-ME)724. 723.723 

724  WRITE  OUTPUT  TAPE  24* 1 1 • J»B < 1 . J ) »B ( 2* J ) »B ( 3* J) »Z .T *B ( 6 * J ) 
GO  TO  702 

723  WRITE  OUTPUT  TAPE  24 • 9» J »B I  1  * J ) »B( 2 * J ) »B< 3* J ) »Z.T *B< 6 * J ) 

702  CONTINUE 

703  IFISENSE  LIGHT  1)704.705 

704  IF(IP) 64 1.54 1*641 

705  S(1)«0.0 

DO  706  J-2.LASTP 
W  «  8  (  3  .  J ) 

DX«B(1»J)-B(1»J-1) 
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DY«8<2*J)-B<2*J-1I 

Si J)*S< J-1)+S0RTF(0X*0X+0Y*DY) 

706  FSU)«B<2#J>*W*B<4»J)*<1.0-W*W)**G1 
IF<IP)708.707.708 

707  SUM-0.0 
GO  TO  709 

708  SUM*FINAL<6»NN) 

709  00  710  J-IGO.LASTP 
K  *  J 

CALL  PARAB<S(K).FS<K).D<1>) 

AOD»POLNT (D<l)tS(K)»S(K+l)) 

SUM  *  ADD  ♦  SUM 
DEL  *  XMASS  -  SUM 
I F (DEL )712»711*710 

710  CONTINUE 
IF(IP)641#54l»641 

711  Fl-1.0 
F2  «  0,0 
GO  TO  717 

712  AREA*ADD+DEL 
K*0 

XLO*S( J) 

XHl*S( J+l ) 

713  K-K+l 

XX ■ ( XLO+ XHI ) *0*  5 
YY*POLNT  <D(1)*S(J)»XX) 

test*area-yy 

IF(TEST)714»716»715 

714  XHI-XX 
IF(K-20)713t716.716 

715  XLO-XX 

IFIK-20) 713»716»716 

716  DX*S( J+l  l-S(  J  J 
F1«(XX-S< J))/DX 
F2«(S<J+1)-XX)/DX 

717  DO  718  1*1 » 5 

718  B<  I*  J*-1)«B!  ItJ)*F2+B<I»  J+1)*F1 
B(6* J+l )  *  1.0/8<4tJ+l) 

DO  719  K*1 *6 

719  WALL!K#ME)*B<K» J+l ) 

LAST* J+l 

IF( IP)720#538»720 

720  IF(N-ME)640*600t538 

800  GO  TO  17 

804  END  FILE  24 

807  STOP  05314 
END 

FUNCTION  POLNT(DfYtX) 

DIMENSION  D<3) 

D2  «  D(2)*0.5 

D3  «  D(3)*0. 33333333 

POLNT  -  (X-Y)*(D(1I-MX+YJ*(D2+D3*X)+D3*Y*Y) 

RETURN 

END 
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SQRTF 

ASINF 

COSRF 

ATANF 

SINRF 

FMA 


FMV 


PARAB 


POL  NT 


OFELD 


APPENDIX  II 
SUBROUTINES  USED 


calculates  square  root  of  the  argument 
calculates  angle  of  argument  (in  radians) 
cosine,  in  radians 
gives  arctan,  in  radians 
sine,  in  radians 

calculates  Mach  number  as  a  function  of  area  ratio, 

(A/A*  or  A*/A),  and  gamma, 
e.g.,  XM  =  FMA  (RR,  G) 
where  RR  =  area  ratio 
G  =  gamma 

calculates  Mach  number  as  a  function  of  Prandtl- Meyer 
angle  and  gamma, 
e.g.,  XM  =  FMV  (R,  G) 
where  R  =  Prandtl- Meyer  angle 
G  =  gamma 

calculates  the  coefficients,  Al,  A2,  and  A3,  of  a  parabolic 
fit  through  three  points, 

YY ( 1-3)  =  FUNCTION  of  XX (1-3) 

i.e. ,  Y  =  Al  +  A2X  +  A3X2 

integrates  the  area  under  the  parabola  Al  +  A2X  +  A3X2. 
i.  e.,  /  ‘(Al  +  A2X  +  A3X2)  dX 
e.g-»  sum  =  POLNT  (A(L),  Xj,  X2) 

calculates  X,  Y,  V/Vmax,  Mach  angle,  flow  angle,  and 
entropy  at  the  intersection  of  a  left  and  a  right  -  running 
characteristic, 

e.g.,  CALL  OFELD(A(l,  i),  B(l,  i),  B(l,  i+1),  G) 
where  A(l,  i)  =  X-coordinate  of  right- running  characteristic 
B(l,  i)  =  X-coordinate  of  left- running  characteristic 
B(L,  i+1)  =  X-coordinate  of  intersection 
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FUNCTION  POLNT ( D*Y  »X ) 

DIMENSION  D( ? ) 

D2  «  D<2)*0.5 

D3  *  D<3)*0. 33333333 

POLNT  ■  IX-Y)*(D(1)+(X+Y)*(D2+D3*X)+D3*Y*Y) 

RE1U9N 

END 
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FUNCTION  FHV 

FUNCTION  FMV  RETURNS  THE  MACH  NO.  AS  A  FUNCTION  OF  GAMMA  P.M.  ANG 

FUNCTION  FMV ( R  *  G ) 

A»SORTF ( < 641.0 J/(G-1.0) | 

B-1.0/A 
Pl-4. 898979 

Rl»A*  ATANF(B*P1)-  ATANFCPl ) 

IF (R-Rl )  1 . 10.2 

1  P2»2. 828427 
GO  TO  3 

2  P2»8. 944272 

3  R2-A*  ATANFCB*P2)“  ATANFCP2 ) 

IF (R2“R)4*11»4 

4  P3»P14(R-R1)*(P1-P2I/(R1-R2) 

IFCP3-1.0) 14* 15 *15 

14  P3-1.0 

15  IF( ABSF (P2“P3 I -0.0005) 12*12*5 

5  R3»A*  ATANF(B*P3I"  ATANFCP3 ) 

Pl«P24(R-R2l*(P2“P3)/iR2-R3» 

IF(P1-1. 0)16.17.17 

16  Pl-1.0 

17  IF(ABSF(P3-Pl)-0. 0005)10. 10. 6 

6  Rl»A*  ATANF(B*P1)“  ATANFCPl) 

P2«P34(R-R3)*CP3“Pl)/(R3“Rl) 

IFCP2 “1.0)18. 19.19 

18  P2-1.0 

19  IFCABSF(Pl-P2)-0. 0005)11. 11. 3 

10  FMV-SQRTFCPl#*24l,0) 

RETURN 

11  FMV-SORTFCP2**24l.O) 

RETURN 

12  FMV«SQRTF(P3**24l.O) 

RETURN 

END 
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C  SUBROUTINE  PARAB 

SUBROUTINE  PARAB! XX.YY.A I 

SUBROUTINE  PARAB  RETURNS  THE  COEFFICIENTS  A  OF  A  PARABOLIC  FIT 
THROUGH  THREE  POINTS  YYI1-3)  «  FUNCTION  OF  XX<l-3) 

DIMENSION  XX<3).YY<31 *A<3> 

XI  >  XX(1) 

X2  »  XX( 2 ) 

X3  *  XX( 3 ) 

FI  «  XI  -  X2 
F2  •  X2  -  X3 
F3  •  X3  -  XI 
D  «  -1.0/(F1*F2*F35 
Y1  *  YY( 1 ) *F2 
Y2  »  YY< 2 1 *F3 
Y3  •  YYI 3 ) *F1 

AU)  «  D*!Y1*X2*X3  ♦  Y2*X1*X3  ♦  Y3*X1*X2> 

A(2 )  »  -D*(Y1*<X2+X3)  ♦  Y2*CX1+X3)  ♦  Y3*<X1+X2)) 

A( 3 )  «  D*< Y1+Y2+Y3) 

RETURN 
END 
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FUNCTION  FMA(RR.G) 

FUNCTION  FMA  RETURNS  THE  MACH  NO.  AS  A  FUNCTION  OF  AREA  RATIO 
RR*, AREA  RATIO  (  A/A*  OR  A*/A  )  ANO  G-GAMMA 

IF (RR  -  1.0)20.21,22 

20  R  «  RR 
GO  TO  23 

21  FMA  -  1.0 
RETURN 

22  R  •  1.0/RR 

23  CONTINUE 

I F ( G-GAMMA ) 1 . 2 . 1 

1  GAMMA »G 

C1«(G*1.0)/(2.0*(G-1.0)) 

C2*( (G*1.0)#0,5)*»C1 
C3*0.5*(G-1.0) 

2  Pl-5.0 

Rl*C2*Pl/<1.0+C3*Pl»Pl)«»Cl 

IF(R1-R)3*10»4 

3  P2-3.0 
GO  TO  3 

4  P2-9.0 

5  R2-C2*P2/<1.0+C3*P2*P2)**C1 
IF(R2-R)6.11,6 

6  P3*Pl+IR-Rl)*«Pl-P2)/(Rl-R2) 

IFIP3-1. 0)14. 15.15 

14  P3-1.0 

15  IF( A8SF { P2-P3 ) -0.0005) 12.12.7 

7  R3»C2*P3/?1.0^C3*P3*P3)**C1 
P1«P2-MR-R2)*(P2-P3)/<R2-R3) 

IF(P1-1.0)16.17.17 

16  Pl-1.0 

17  IF  ( ABSF(P3**P1 ) -0.0005)  10.10 ,8 

8  Rl*C2*Pl/< 1.0+C3*P1*P1)**C1 
P2«P3*<R-R3)*<P3“Pl)/tR3-Rl ) 

IF (P2-1.0) 18.19.19 

18  P2-1.0 

19  IF(ABSF(Pl-P2)-0. 0005)11. 11. 5 

10  FMA-P1 
RETURN 

11  FMA»P2 
RETURN 

12  FMA*P3 
RETURN 
END 
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SUBROUTINE  OFELD  CA.B.C.G) 
DIMENSION  A C  5 ) »B<  51 »C<  5 ) 

XI  ■  A<1) 

Y1  «  A  f  2 ) 

W1  *  A 1 3  > 

SA1  ■  AC  A) 

T 1  *  AC  5  ) 

W13  «  W1 

X2  »  BCD 
Y2  «  BC2) 

W2  *  BC 3 ) 

SA2  *  BC A) 

T 2  r  BC5I 
W23  «  W2 

I F ( SENSE  SWITCH  A>91*92 
CONTINUE 

WRITE  TAPE  12  #OASH 
WRITE  TAPE  12»CACJ)*J*1»5) 
WRITE  TAPE  12t«BC J) 

WRITE  TAPE  12»BLANK 

CONTINUE 

DT  «  T1  -  T2 

OW  «  W1  -  W2 

C3  «  0,0 

I  «  -1 

STl  ■  SlNRFCTlf 
HST1*  STl 
CT1  »  C0SRFCT1) 

HSA1«  SA1 

CA1  ■  SORTFC1.0  -  SA1*SA1» 
HCA1*  CA1 
ST2  «  SINRFCT2) 

HST2«  ST2 
CT2  *  COSRF  C  T 2  1 
HSA2*  SA2 

CA2  «  SORTFC 1*0  -  SA2*SA2I 
HCA2*  CA2 

SINA  »  ST1*CA1  -  CT1*SA1 
SA  -  SINA 

COSA  «  CT1*CA1  ♦  ST1*SA1 
CA  ■  COSA 

SINB  *  ST2*CA2  ♦  CT2*SA2 
SB  ■  SINB 

CO  SB  *  CT2KA2  -  ST2*SA2 
CB  *  COSB 
60  TO  AO 

ST3  «  SINRFCT3) 

CT3  «  C0SRFCT3) 

CA3  ■  SORTFC 1.0  -  SA3*SA3) 
HI  •  ST3*CA3 
H2  «  CT3*SA3 
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H3 

m 

CT3*CA3 

H4 

m 

ST3*SA3 

SINA 

w 

(SA  ♦ 

H1-H2  >*0«  5 

COSA 

* 

(CA  ♦ 

H3+H4)#0,5 

SINB 

w 

(SB  ♦ 

H16H2>*0.5 

COSB 

m 

(CB  ♦ 

H3-H4)*0.5 

ST1 

9 

(HST1 

♦  ST3>*0.5 

SA1 

m 

( HSA1 

♦  SA3)*0«5 

CA1 

9 

(HCA1 

♦  CA3>*0.5 

ST  2 

9 

(HST2 

♦  ST3)*0*5 

SA2 

9 

(HSA2 

♦  SA3 ) *0.5 

CA2 

9 

(HCA2 

♦  CA3)*0.5 

W13 

9 

(W1  ♦ 

W3)*0.5 

W23 

9 

(W2  ♦ 

W3>*0.5 

C3 

9 

W3 

CONTINUE 

E  ■  SINA*COSB 

F  «  COSA*SIHB 

0  «  F  -  E 

X3  «  <F*X2  -  E*X1  ♦  COSA*COSB*( Y1-Y2) )/0 
Y3  «  (F*Y1  -  E*Y2  ♦  SINA*SI NB*( X2-X1 ) ) /D 

OX  «  X3  -  XI 

OY  «  Y3  -  Y1 

IFIABSF(OX)  -  ABSFIDY) >41.42.42 

41  PI  «  OY/SINA 
60  TO  43 

42  PI  «  DX/COSA 

43  OX  «  X3  -  X2 

OY  «  Y3  -  Y2 

IF(ABSFIDX)  -  ABSF(DY) >44.45.45 

44  P2  «  OY/SINB 
60  TO  45 

45  P2  «  OX/COSB 

Pi  «  STl*SAl*Pl*2.0/{ Y1+Y3) 

P 2  «  ST2*SA2*P2*2.U/<Y2+Y3) 

50  Ql  *  CA1/(SA1*W13) 

02  ■  CA2/(SA2*W23) 

Q  «  02  ♦  Ql 

T 3  «  (Q2*(P1  ♦  T1  +  Q1*DW)  ♦  Q1*(T2-P2 ) ) /Q 

W3  «  (P1+P2  ♦  Q1*W1  +  Q2*W2  +  DT)/Q 
WW  *  W3*W3 

SA3  «  SQRTFU1.0  -  WW)/(6*WW)  ) 

I  «  I  ♦  1 

IFJSENSE  SWITCH  4)93*94 

93  CONTINUE 

WRITE  Tape  12.P1.P2.Q1.02 
WRITE  TAPE  12.X3.Y3.W3.SA3.T3 

94  CONTINUE 

IFC I >52.51.52 

51  TOLD  ■  T3 

T3  «  ( T 3  ♦  < T 1+T2 )*0*5 ) *0*5 
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60  TO  30 

52  IFtABSFt T3-T0LD)  -  0.00001 )60. 60*61 

60  IF(ABSFIC3-W3)  -  0.00001170*70*61 

• 

61  I F  C 1-40)62  *63*62 

62  TEMP  *  T3 

T3  *  <T3  ♦  T0L0)*0.5 
TOLD  *  TEMP 
60  TO  30 
C 

63  SENSE  LI6HT  1 

70  CC1)  *  X3 

C  C  2  >  *  Y3 

CCS)  «  W3 

C<4)  •  SA3 

CC5)  •  T3 

80  RETURN 

DASH  »  DASH 
BLANK  ■  BLANK 
END 
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APPENDIX  III 

FIELD  POINT  CALCULATION 


The  partial  differential  equations  for  supersonic,  axially  symmetric, 
irrotational  flow  are  reduced  to  the  following  ordinary  differential 
equations  by  the  method  of  characteristics  (Ref.  1): 


1.  = 
dX 

tan  (0 

-  a) 

(right-running) 

2.  -43L  = 

dX 

tan  (0 

+  a) 

( left-running) 

3.  id  + 

cot  a 

(df)  - 

sin  0  sin  a  dY 

W 

sin  (  0  +  a  )  Y 

4.  id  - 

cot  a 

(df)  + 

sin  0  sin  a  dY 

V 

sin  (  0  +  a )  Y 

0  (right-running) 

0  ( left-running  ) 


Equations  (1-4)  are  solved  by  finite  difference  methods  in  the 
following  form: 


Yc  *  sin  ( 0  -  a ) 

Xc-X.  cos  (  0  —  a  ) 

Yc  ~Yb  s in  ( 0  +  a) 

Xc  -X.  cos  ( 0  +  a) 

1.  Xc  •  »i»  (0  -  a)  -  Yc  •  cos  (0  -  a)  =  XB  •  sin  (0  -  a)  -  Y,  •  cos  (0  -  a) 

2.  Xc  •  sin  (0  +  a)  —  Yc  •  cos  (0  +  a)  =  Xb  •  sin  (0  +  a)  -  Yb  •  cos  (0  +  a) 

Solving  the  system  Eqs.  (1)  and  (2),  we  have 

X,  •  sin  (0  —  a)  —  Y,  •  cos  (0  —  a)  —  cos  (0  —  a) 

X|,  •  sin  (0  +  a)  —  Yb  •  cos  (0  +  a)  —  cos  ( 0  +  a) 

•  sin  (0  —  a)  —  cos  (0  —  a)  I 

sin  (0  +  a)  —  cos  (0  +  a)  | 

sin  (0  —  a)  X,  •  sin  (0  —  a)  —  Y,  •  cos  (0  —  a) 

sin  (0  +  a)  Xb  •  sin  (0  +  a)  —  Yb  •  cos  (0  +  c) 

|  sin  (0  -  a)  —  cos  (0  —  a) 

|sin(0  +  a)  —  cos  (0  +  a) 


1. 

2. 


3.  0C  -  0„  + 


cot  a 

W 


(Wc  -  Wa)  - 


sin  0  sin  a 
sin  (  0  -  a  ) 


il 

Y 


0 


4.  0C  -  0b  -  (Wc  -  Wb)  + 


sin  0  sin  a  dY 
sin  (  0  +  a  )  Y 


0 
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APPENDIX  IV 
DATA  PRESENTATION 


This  appendix  given  a  sample  of  the  output  obtained  by  the  methods 
previously  outlined.  For  each  of  the  two  principle  sections,  there  is 
given  a  page  of  wall  points,  a  page  of  axis  points,  and  several  pages 
of  characteristic  lines  including  the  first  and  last  characteristics  in 
the  respective  regions. 

As  indicated  in  Section  2.  3,  each  characteristic  originating  on 
line  EF  has  an  initial  mass  integration  value  which  is  computed  at  the 
first  point  of  that  characteristic.  This  initial  value  of  mass  appears  in 
the  output  as  the  first  value  in  column  seven. 
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SAMPLE  CALCULATIONS  TFRCAT  REGION 


INPUT 

NUMBER  CF  PCINTS  CN  FIRST  CHARACTERISTIC 
NUMBER  CF  PCINTS  CN  AXIS  (N)  =  20 
GAMMA=  1.4000 

INFLECTION  ANGLE  (ETA)*  12.000  DEGREES 
COORDINATE  CF  INFLECTION  POINT  (XC>= 
FACTORS  IN  V/VMAX  TC  CETERMINE  X  t  X 1 )  = 

{ X2 )  = 

POWER  GOVERNING  AXIS  DISTRIBUTION  (P)=  2 
COEFFICIENTS  IN  TERMS  OF  V/VMAX 
C(l)=  4.0824825-01 
C(2  )*  5.5127312-01 
C (3 )=  1 .0438091-02 
C<4)*-3. 6969777-01 
C ( 5 ) *  2.2921343-01 
C(6)*-4. 2703583-02 


(M)*  15 


OCOOO 

.47270 

1.91607 

.00000 
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